International Journal of Advanced Engineering Research and Science (IJAERS) 
https://dx.doi.ora/10.22161/iiaers.5.11.8 


[Vol-S, Issue-11, Nov-2018] 
ISSN: 2349-649S(P) / 2456-1908(0) 


Remote Sensing Satellites Planning System 

Amr Qamar 1 , Essam Eldin S. 2 , Khaled M. Badran 3 

'Eng., MID, Egypt 
2 Dr., MID, Egypt 
3 Dr., MTC, Egypt 


Abstract — A Remote Sensing Satellites Planning system 
(RSSP)for satellite constellations is responsible for managing 
these satellites by assigning the imaging tasks to each satellite 
in the constellation such that the loads are balanced and the 
resources are well used. The proposed system can be used 
with heterogeneous constellations that consist of satellites 
whose different specifications, different orbits' types and/or 
different payload types. This problem is a combinatorial 
optimization NP-hard problem modeled in this paper as a 
Constraint Satisfaction Problem using the Constraint 
Programming Technique. The output plan is obtained using 
one of three objective functions (gain maximization, area 
maximization, and image quality maximization) using four 
search algorithms (simulated annealing, hill climbing, tabu 
search and late acceptance) and different planning horizons 
(one track, one day and one month). 

Keywords — Agile satellites, remote sensing, planning, 
constraint satisfaction problem, heterogeneous constellation. 

I. INTRODUCTION 

The output plans from a RSSP system should satisfy the 
customers by performing the maximum number of their 
requests in addition to optimize the usage of the constellation's 
resources taking into consideration all the system and 
operational constraints to assure obtaining a feasible output 
plan. 

The paper introduces the related work in Section(II). Then 
the proposed RSSP system formulation in Section(III) is 
discussed. The system architecture is introduced in 
Section(IV). Section(V) clarify the system performance and 
introduces some comparisons with other systems from CPU 
time point of view. The final Conclusion and 
recommendations for future work are introduced in 
Section(VI). 

II. RELATED WORK 

The constraint programming technique were used by two 
references, in 2015 Erik Demeulemeester et al [1] used non- 
agile earth observing satellites in their system The used 
planning horizons were multiple tracks and one day. They 


used branch and price search algorithm and column generation 
heuristics. While in 2002 M., Verfaillie, G., Jouhaud, F., 
Lachiver, J., & Bataille and N. Lemaitre [2] used single agile 
satellite with the Cieedy search algorithm and the simulated 
annealing. They describe their system for the illuminated half 
of the track. 

III. PROBLEM FORMULATION 

The problem model can be partitioned into the following four 
parts: 

1) Input Data 

For the set of input requests R. there exists r £ R, whose gain 
G- and surface area A r . Let I be the set of images obtained 
from R by the geometric cutting up process. For each image i 
£ I: 

Ei : earliest shooting time, Li : latest shooting time, Di : 
duration of shooting, Ai : surface area. For each possible pair 
of images (i, j), Mij is defined as the transition time between 
shooting the two images consecutively. B is defined as the set 
of pairs of images (i, j) such that i and j are images for the 
same strip with opposite pitch angles and S be the set of pairs 
of stereoscopic images. 

2) Decision Variables 

a) Xi ... is One if the image i £ I is selected, and Zero 
otherwise. 

b) ti ... the shooting start time of image i if selected. 

3) Constraints 

The following constraints have been implemented in the 
proposed system: 

a) Observation Time Window Constraint: Vi £ I: 


(Xi = 1) => (Ei < tj < LO (1) 

b) Transition Time Constraint: V (i, j) £ I : (X i; Xj = 

1) => (tj + Dj + Mi, < tj) (2) 

c) Request End Time Constraint: Vi £ I, I c r,r £ 

R: (X; = 1) => (tj < end r ) (3) 

d) Mono Image is Shot Once Constraint: V (i, j) £ B : 

Xi + Xj < 1 (4) 

e) Stereo Image Constraint: V (i, j) £ S : Xj = Xj 

(5) 
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f) The Stereo Pair Pitch Angles Constraint: V (i, j) £ 

S ; Pitch) = — Pitch j (6) 

g) The Stereo Pair Roll Angles Constraint: V (i, j) £ S : 

Roll) = Rollj (7) 

h) The Stereo Pair Satellite Constraint:V (i, j) £ S : 

Satellite) = Satellite^ (8) 

i) Onboard Memory Constraint: V planning horizon, 

Vi £ I: E)X) < N max (9) 

Where Nmax is the maximum allowable number of 
selected images for this planning horizon. 

j) Payload Duty Cycle Constraint: V (i,j) £ I=tj — 

(t)+ D) + M),) > t P (10) 

Where tp is the payload duty cycle or the so called 
technological break. 


[Vol-S, Issue-11, Nov-2018] 
ISSN: 2349-649S(P) 12456-1908(0) 

4) Objective Functions 

The three objective functions implemented in the proposed 
systemare defined as follows: 

a) Gain Maximization function: V i £ I : f(x) = 

max(Z)Gi*Xi) (11) 

b) Area Maximization function: V i £ I : f(x) = 

max(Z)A) * X)) (12) 

Where A) = D) * 

W s (13) andW s ...is satellite's swath. 

c) Image Quality Maximization function: V i £ I : 

f(x) = max (Z)(R worst - R)) * X)) (14) 

R)... is the resolution of image i £ I and R worst ... 
resolution at maximum roll and pitch angles. 


IV. PROPOSED SYSTEM ARCHITECTURE 


The operation through the proposed system architecture in Fig. 1 is explained as follows: 



Requests database will contain the input requests with their 
data. The Geometric cutting up process will produce images in 
Images database that is fed to the Orbit Propagator with the 
satellites in Satellites database and the selected planning 
horizon. The Orbit Propagator outputs the imaging 


opportunities for each image. The Opportunities database will 
contain all the imaging opportunities for all the images with all 
the satellites. The operation in the RSSP system begins with 
modeling the input via the Modeler to be sent to the SOLVER 
that is configured using a Solver Configurator. To build the 
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Constraints database and the Objectives database, the 
Constraint Creator and the Objective Builder is used 
respectively. The Score Comparator takes the created 
constraints and the built objectives as inputs. The SOLVER 
searches for solutions in the search space and selects feasible 
solutions during the algorithm-running lifetime. Each time, the 
selected solution is score calculated via the Score Calculator. 

This score is compared with the previous scores via the Score 
Comparator until reaching the optimum score. 



V. TESTS AND RESULTS ANALYSIS 

To determine the system’s performance, it is tested for three 
case studies with different problem sizes and measure the CPU 
time (efficiency) and the Normalized Score (quality). The first 
case study consists of (1) satellite, (1) track planning horizon 
and (25) targets while the second is (2) satellites, (1) day and 
(120) targets and the third one uses (3) satellites, (1) month 
planning horizon and (2500) targets. Fig. 2 illustrates the 
system’s behavior in the (3) objectives using the (4) search 
algorithms and measures the (2) metrics. 


It is clear from the figure that almost all the results lie in the 
left bottom part of the graph which means that these results 
has high normalized score and low CPU time. This 
emphasizes that the system is qualified and efficient. We 
compared our results, from the CPU time point of view, with 


those published for similar systems according to the problem 
size. 

The comparison introduced in TABLE I clarifies that the 
proposed system performs very well compared with the listed 
references. 
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Table I 

COMPARISON WITH OTHER SYSTEMS 


in 

Sat's 

no. 

The Reference Results 

Proposed SystemResults 

-u. 2 

D. 

Ref. 

# of Targets/ Imaging 
Opportunity 

CPU Time 
(seconds) 

# of Targets/ 

Imaging Opportunity 

CPU Time 

(seconds) 

i 

i 

[3] 

30 

153.2 

25 

32 (Max) 

2 

2, 3 

[1] 

100 (3 Satellites) 

486.12 

120 (2 Satellites) 

304 (Average) 

3 

3 

[4] 

250 Imaging 
Opportunities 

2995.5 

86195 Imaging Opportunities 
(2500 Targets) 

2532 (Max) 


VL CONCLUSION AND FUTURE WORK 
RECOMMENDATIONS 

The proposed system is implemented and tested to be used 
with different target types, different planning horizons and 
different constellation types. The system is designed with 
many constraints and four search algorithms. Three different 
objectives for creating the output plan are used. The tests 
results are relatively good compared with some other similar 
systems. 

It is recommended for the future work to include more 
objectives and use other search algorithms. 
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